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Abstract 

We discuss the foundations and extend the range of applicabihty of the widely used Kubo- 
Greenwood formula (KGF) for the electronic conductivity. The conductivity is derived from the 
current density, and only the probability amplitude rather than the transition probability is used. It 
is shown that the contribution to the conductivity from degenerate states in a low or zero frequency 
external electric field and the contribution from states near resonance with a finite frequency 
external field are finite. The improved conductivity expression does not include the familiar "energy 
conserving" delta function, and no artificial broadening parameter for delta function is required 
for the DC limit. We explored two methods of computing current density. We discuss the role of 
many-electron statistics in computing the conductivity in single-particle approximations, and we 
show that the conventional KGF is due to the contribution from single-particle excited states. 

PACS numbers: 71.15.Pd, 72.10.Bg, 72.20.-1 



I. INTRODUCTION 



The Kubo-Greenwood formula (KGF) has been widely used with great success to calculate 
the electrical conductivity from first principles simulationsi'^i'^i^i^'^. Despite the universal use 
of KGF, it is worth pointing out that it has certain limitations, and for some applications, 
improvements are possible. The aim of this paper is both to derive the KGF and more 
fundamental transport formulae from first principles, and to point out cases in which use of 
the KGF can be problematic, with a full explanation of the origins of the difficulty. New 
formulae are presented which circumvent some of these difficult cases. 

Greenwood's derivation of the conductivity used the transition probability between two 
single-electron states in an oscillating external field, in such a way that the interaction time 
must be long enough to assure that the transition probability is well-defined. On the other 
hand, to make perturbation theory applicable, the interaction time should be short^i^i^ii^. 
For a large system, in which the energy spectrum is continuous, these two conditions are in 
conflict. The usual time-dependent perturbation theory cannot be applied in two cases: (1) 
degenerate states and (2) if two groups of states are in resonance with an oscillating external 
field. Both circumstances are common in a macroscopic system. In addition, the energy 
conserving 6 functions in KGF require an artificial broadening parameter when implemented 
numerically. 

To avoid these difficulties, in Section II, both direct current (dc) conductivity and alter- 
nating current (ac) conductivity are derived from a new expression for the current density. 
Only the probability amplitude (not the probability itself) enters, and the 5-function in the 
KGF disappears in the improved expression of conductivity. The new expression may be 
reduced to the well-known result for the conductivity by applying the Boltzmann equation 
with the relaxation-time approximation for a crystal. For a static field or oscillating field, 
the zeroth-order contribution from a group of degenerate states is shown to vanish, and the 
first-order contribution of degenerate states is finite (Appendix A). We also show that the 
zeroth-order contribution from two groups of states which are near or in resonance with an 
oscillating external field is zero. The first-order contribution of the two groups of resonant 
states is finite. (Appendix B). By invoking the single-particle approximation at different 
stages of the derivation, one is led to slightly different results. In section III, we use many- 
body perturbation theory to compute the conductivity for an intrinsic semiconductor, metal 
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and doped semiconductor. It is shown that at T=OK, the dc conductivity of an intrinsic 
semiconductor is zero, a well-known consequence of many-electron statistics. 



II. THE CURRENT DENSITY 

If the motions of nuclei are treated classically (as in most ab initio MD codes), the average 
energy H of the electron + nuclei system in an electromagnetic field described by vector 
potential A and scalar potential </? is given by 

^ = J d^NA{w^-WM}(^i ■ ■ ■ riVe)^e^{Wi...w^}(ri • • • fivj (1) 

a a a,/3(>o) 

where g^, M^, Pa and are the effective charge, mass, canonical momentum and position 
vector of the a^^ nucleus. V^Wa, W/3) is the interaction between the a^^ nucleus and the 
nucleus. 

H'e = E[^(P^- - ^A(r,))^ + e(^(r,)] + ^ V{rj, r,) + V{r„ (2) 

i=l jk ja 

is the electronic Hamiltonian in the external electromagnetic field. V{rj, r^) is the interaction 
between an electron at and another electron at r^, V{rj,Wa) is the interaction between 
an eletron at and the a^^ nucleus at W^. The motion of the electrons is determined by 

A{Wi---Wjv-}'^'"i ■ ■ '^A^e) is the many-electron wave function of H'^ for a given nuclear config- 
uration {Wi ■ ■ ■ Wjv") belonging to eigenvalue E'f . We use He to denote H'^ when 
external field does not appear, A{Wi - WAr} the eigenfunction of Hg belonging to eigenvalue 
■^{w wj^} ■ Hereafter we use a symbol with prime to denote a quantity when external field 
appears, the corresponding symbol without prime to denote the quantity in zero field. The 
velocity of the 7**^ nucleus is determined by 

where V^, = P^^^'^/M^, is the velocity of the a^^ nucleus, P^^^'^ = Pa - gaA(Wa) is the 
mechanical momentum of the a^^ nucleus. After some manipulations, P-y = 



as 

M,V, = g,[E(W,) + V, X B(W,)) - J] ^^^^^^ - (5) 

In the MD formulation, the positions of the nuclei are functions of 'time' (MD steps). 
The initial positions of nuclei are given from an initial configuration, the initial velocities of 
nuclei are assigned in some way. The electronic wave function A{Wi - WAf}(j"ii^2r3 ■ • ■ r^rj is 
calculated from the configuration {Wi ■ ■ • W^v"}? the forces on each nucleus is then calculated 
from A{Wi- -Waa}(i'iI"2I'3 ■ ■ ■ fAfe)- The position and velocity of a nucleus in next step are 
calculated from the length of the time-step, acceleration and the velocity in last step^. 

According to the principle of virtual work, for a given state A', the microscopic electric 
current density at point r 

= iVe— / rfr2rfr3---t/r^,(A'Vr,A'*-A'*Vr,A')--A(r)n|w,...w^}(r) + $^gaV,5(r-W„) 
where 



'^{Wi 



.,...w^}(r, t) = NeJ drsdrs ■ ■ ■ drNA{-w^-w_^}irr2r-i ■ ■ ■ rN,)^[w^■■■w_^}i^r2rs ■■■ynJ 

(7) 

is the number density of electrons at r for a given nuclear configuration {Wi---W^}. 
Eq.([6]) is the response of the electrons+nuclei system to the external field; the first two 
terms are due to electrons, and the last term is due to nuclei. The measured macroscopic 
current density at point r spatial average of Eq.(l6]) over a region fly. centered at r: 



j(r) = — / rfsj„(s) (8) 

The linear size L of Qr satisfies: a << L << A, where a is a typical bond length, A is the 
wavelength of external field or other macroscopic length scale. Eq.([H]) is the usual current 
density defined for an infinitesimal arear^*^. 

Using the single-electron approximation to separate variables in Eq.([3]), we obtain the 
equation satisfied by the single-electron wave function x'l'- 

KxKt) = Elx'iir), K = ^(P - eA(r))2 + e^{r) + f/(r, {WJ) (9) 
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where /i^ is the single-electron Hamiltonian in an external field, U is the single-electron 
potential due to nuclear configuration {Wq}. ha, and Ei are the corresponding quan- 
tities when external field does not appear. They are the Hamiltonian, eigenfunctions and 
eigenvalues as in density functional theory (DFT), or other single particle theories. 

The current density due to electrons can be computed as following. At finite temper- 
ature T > 0, the system can be in the ground or excited states. The electron current at 
temperature T comes from both the various excited states and the ground state: 



iheN, 



(10) 

where 



e- p 



is a Ai'e— electron state, P is a permutation on objects (riSzi; r2Sz2', rss^s; ■ ■ ■ ; r]\i^SzNe), 
(5p = 1 if P is an even permutation, 6p = —1 if P is odd. Of course h,l2,--- Jn^ are 
distinct. Because any observable like j*^ is bilinear about Aj^^^...^^ , the order of rows and the 
order of columns in A;^;^...;^ do not matter. We only need to maintain a fixed order in all 
intermediate steps of calculation. The sum is over all possible choices of A^e single-electron 
states. The arguments of A' are (s, r2,r3,-- - ,rNj, to save space the spin variables are 
abbreviated. 

WU-i..=UU-ijZ', Z'= Yl ^.z....... ^.z...... =exp[-(i?U...,,^-£;^)/(A;^T)] 

(12) 

is the appearing probability of state AJ^^^...;^ . Eq is the energy of A'e— electron ground state. 
When no field is applied on the system, macroscopic current does not appear in any state 
^hh-iNe- '^^^ current density from electrons reads: 

= d^jdr2dv,...dv^^ Y nz....j(AU...^,V.A-,..,^^-A-^..,^^V.A;^,,..,, 

- (Az,z,...z^,VA;z,...z^, - Ar,i,..,^^V.A,,i,...^^J](sr2r3---r^J (13) 

For low temperatures, {xhl ~ eE ■ r\xii) << ksT is not satisfied. Linearizing Wl_^i^,,,i^ 
about field E is not legitimate (cf. Eq. f[T^ ): current density is not necessary linear 
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about field, the dependence of conductivity on field is intrinsic at low temperature for 
semiconductors. If temperature is not too low ((x/J ~ ^Er\xh) << ksT), we may expand 
W'[] in Eq.^ to first order of field 

a a a 

The last equal sign used the obvious fact We=o = 0: no macroscopic current exist when 
external field vanishes. With the help of Eq. (fT4l) . Eq. (fT3!) is simplified to 



J w 



iheNe 



- (A«i/2.../^,VsA;, - Ali^...i VsAhi,-i^J]isr2rs ■ ■ - r^J (15) 



where 



(16) 

are the corresponding quantities without external field. In the single-particle approximation 
W,,,,..,,^=l[f{Ea /(E^J = -^^--^^.-^ (17) 

where n is chemical potential at given temperature and shape of the interested body. 

The current density (ITSl) and the conductivity deduced from it are just for one MD step. 
To include the the thermal vibrations in a material, one must average the conductivity 
over many MD steps. Only the averaged conductivity may be compared to the experimen- 
tal observations where the material changes its configurations with time through thermal 
vibrations. This observation is valid for solids, liquids and molecules. 

The idea of linear response^ can be applied in two different ways: (1) first express AJ^^^...^^ 
with single-electron wave functions x' and effect the multiple integral J dr2dT^ ■ ■ ■ dr^^. Then 
view x'la cLS correction of Xia under perturbation — eE ■ r. (2) view AJ^^^...^^ as correction of 
Ajjij-.-Ziv^ under perturbation — J2m=i sE-r^. Then effect multiple integral / dr2dr3 ■ ■ ■ dr^^. 
Most discussion about the KGF-'^'^ii^ is based on method (1). The role of many-electron 
statistics is displayed more explicitly in method (2). Both schemes express the conductivity 
in terms of single-electron states and corresponding eigenvalues. In the remainder of this 
section, we will use method (1) and compare with previous result. Method (2) will be 
analyzed in the next section. 
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A. DC conductivity 



Applying Eq. (fTT!) and working out / dY2dY^ ■ ■ ■ dr^^, Eq. (fT5!) leads to: 
ihe f 

(18) 

the argument of all single-electron functions is s. 

In a static electric field, the nuclei and the bound electrons are pushed in opposite direc- 
tions. These lead to a static deformation of the material. Since a static electric field does 
not produce any net velocities of nuclei, the 3rd term in Eq.([6]) is zero. A static electric 
field is solely determined by scalar potential v'(r), which means A(r) = 0. The 2nd term in 
Eq-© vanishes. The interaction with an electron at r is 

Hf^ = e^(r) = -eE • r (19) 

At this point, let us assume all the single-electron states in Eq.( |T8i) are non-degenerate. The 
case of degenerate states will be discussed later. From first-order perturbation theory, the 
change x'c^^ fhe single electron wave function due to the external field is 

Xc and Ec are the single electron wave function and the corresponding eigenvalue without 
external field. We should emphasize that voltage is proportional to the distance between 
two points, and so too is the interaction. The change in states cannot be described by the 
perturbation result Eq. fl20p . Except for very weak field, one must use WKB method rather 
than perturbation theory. In this work, let us limit ourselves to very weak field. Substituting 
Eq. fl20l) into Eq. ffTSj) . and only keeping the terms linear with external field, one has 

{{Xi\eE ■ r\xiJ{x:ysXi - Xi'^sXiJ - {Xi\eE ■ r\xiJ*{xiysXi - Xi'^sXiJ} 

The sum over /(t^ l^) is not restrict to (/1/2 ■ ■ -^atJ; it extends to all single particle states. 
By means of the definition of conductivity a^i, 

U = Y (^fjivEy, /X, z/ = X, z (22) 



the dc conductivity is 

E ^W.„.i:E^'™tol^.lxo/Ato^-xL|j) (23) 

In a large system, the matrix element of position operator is not well defined. Making use 
of ^ 

one can change the matrix element of position operator into the matrix element of momentum 
operator-. 

Current use of the KGF is amounts to assuming that beside /„, other single-electron 
states in - ■ - In^ are occupied. Only the factor f{Ei^ is left. Thus the sum over 
various choices of {Z1Z2 ■ ■ - In^ can be ignored if one extends the sum over a to all possible 
single particle states. 

In parallel with Greenwood's work for ac field, Luttinger has derived an expression for 
static field by adiabatically introducing the interaction.— Transition probability rather than 
the amplitude of probability was used. Eg. (123!) does not obviously display a feature of an 
intrinsic semiconductor: dc conductivity vanishes at zero temperature. In addition, due to 
the use of the single particle approximation before applying perturbation theory in Eq. (!23l) . 
one cannot exclude coupling between two occupied states. These faults can be cured in time- 
dependent perturbation theory or by applying perturbation theory directly to the many- 
electron wave function. 

If there is only one group M degenerate single-electron states {xda-, ^ = 2, ■ ■ ■ , M) in 
^hM,- ,iNe^ we first form correct zeroth order wave functions 

xr = E^'^^'^.'^'^." ^,^' = l,2,---,M (25) 

a' 

the secular equation satisfied by Cd^d^, is 

Y^iVd^d^, - edd^dJCd^d^, = 0, Vd^d^, = I dvxU-^^ ■ (26) 

The perturbation matrix {Vd^d^,) is Hermitian, it can be diagonalized by a unitary transfor- 
mation iCd^d^i)- Therefore 

E(^2>xr* - x'ifvxr ) = $:(x..vxi - vx.j (2?) 



According to Eq. (fT8|) . the zeroth-order contribution of the degenerate states to conductivity 
is zero. This is consistent with usual experience: an electron is not accelerated along the 
direction of the field when it transits between states with same energy, and thus makes no 
contribution to the conductivity. 

The first order correction to "x^^ is^ 



where 



Sd. = j rfrxX(-eE • r)x',^, V,d. = j rfrx^-eE ■ r)xf , V^^,, = j rfrxJ°?*(-eE ■ r)xk 

(29) 

k indexes non-degenerate states. The first sum in Eg. (1280 runs over all states which are not 
degenerate with {xd^i = 2, ■ ■ ■ , M). 

Making use of Eq. (1281) , the macroscopic current density Eq. (JT8i) becomes 



' 1.1 ^—1 tl-Lt \ 'Q ' "'Si^/ 



M 



+EEi-%!^%^/ *(x.v.xr - xr-v.x.) (30) 
+E E ""li-^E Jf^i /„ -^Kxrv^x^rNx^^-xr-)} 

The sums over I and k are not restrict to (/1Z2 ■ ■ ■ ^Ne)- The first term in bracket is the con- 
tribution from coupling among non-degenerate states, the third term is contribution from 
coupling among M-fold degenerate states, the second term is contribution from coupling be- 
tween member of non-degenerate states and the M-fold degenerate states. All denominators 
are non-zero. 

Because Xd^\ and Vkd^ are functions of field, the 2^^*^ term and the S'''^ term in Eq. (l30|) 
do not exhibit a simple linear relation with field: each contribution to conductivity is field- 
dependent. The generalization to the situation in which there are several groups of degen- 
erate states in A^^^ij,--- ,'iVe straightforward. 

The massive degeneracies in the N^-electron states Ki^^i^^... ^ij^^ do not cause any trouble: 
they are counted by the sum over all possible ways of picking up single-electron states. In 



a crystal, the degeneracies in the single- electron states may be caused by the high symmetry 
of the lattice. Sometimes accidental degeneracies (band crossing along a symmetry axis or 
a symmetry plane not compelled by symmetry) also occur^"^), as well as additional degen- 
eracies produced by time reversal symmetryi^ii^. In a disordered system (liquid, amorphous 
solid and some molecules), only the degeneracies produced by time reversal symmetry are 
left. From the perspective of computing the conductivity, a disordered system has less degen- 
eracy, and is simpler to treat than a crystal. Comparing the huge number of non-degenerate 
states (most of original degeneracies in the unit cell are removed by the interaction with 
first, second and third coordination shells), the number of degenerate states in each manifold 
is quite small. The third term in Eq. fl30l) . the coupling among degenerate states, is small 
compared with the first two terms. If we neglect the field-dependence of x'd^^ in the second 
term of Eq. fl30l) and use the original Xd^, one does not need to distinguish degenerate states 
and non-degenerate states. Eq. (!23!) could be used to compute conductivity by neglecting 
the coupling among degenerate states. 

B. AC conductivity 

The macroscopic current density in an oscillating field is 

•fc P 

«---*'=2S^X^^ 5: E (31) 

In an ac electric field E = Eq cos out, the interaction of an electron at r with field is 

HfUt) = Fe-'''^ + Fe**", F = -^er ■ Eq (32) 

Since the region Qj. (employed to compute the spatial average) is much smaller than the 
wavelength of the field, the position dependence of field is ignored in Eq. fl32l) . The wave 
function x'd't) ^in external field can be computed from time-dependent perturbation theory 

X'M = Xc{t) + J2 «^(^)Xde-^*^'^/^ Xc{t) = e-*^=/\,(r) (33) 
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where ad(t) satisfies 



dt 

where 



F^^ = J dvxlFxc, ^dc, = \{Ed - E,,) (35) 
We assume initially only state Xc is occupied and other states are empty: a^^ (t = — oo) = 6cci- 



1. non- degenerate states 



If all the states are not degenerate, Eq. (lMI) is simplified to 



dt h "'^ h 

The solution of Eg. (1361) is simply a time integral 



flKUJ^c — UJ — lO) h[UJdc + UJ — lO) 

The change in wave function Xc due to the ac field is 

^}'\t) = 5^ (1 - n,)a,(t)x.e-^*^^/'^ (38) 

where is the occupation probability of single-electron state Xd- The "un-occupation" 
factor (1 — nd) is implicitly assumed when we simplify Eq. llMl) to Eq. (l36|) : initially state 
Xd must be empty. In an intrinsic semiconductor at zero temperature, the valence band is 
fully occupied and conduction band is empty. Eq. (1551) indicates that only conduction states 
can couple with the valence states. This implies that the dc conductivity of an intrinsic 
semiconductor at zero temperature is zero. 

With the help of Eq.([33D, Eq.([37D and Eq.([38D, Eq.([3l]) becomes 

{{Xd\AxiS [-( V + 7 iXicy^xk- VsXz. ) (39) 

+ {Xd\r\xiJ[^^ T + T- — T](XdVsXL - xiysXd)} 
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Separate out cos cut terms in Eg. (1391) : they are in phase with the external field. The real 
part of conductivity then reads: 



J ^ /J- 

When = Eq. lHOl) is reduces to Eq. (l23ll except for the factor (1 — n^), as it should be. 

To compare with the result from Boltzmann equation, let us consider a crystalline metal. 
The semi-classical current density is given by^S 

^M^)9i^) (41) 

where the integral is over the first Brillouin zone, (^(k) is the non-equilibrium distribution 
function under external field, in relaxation time approximation^ 

^7(k) = /(k) + |^eE-vr[E(k)] (42) 

where r[i?(k)] is the energy-dependent relaxation time, /(k) is the Fermi distribution func- 
tion f[E{k)\. The conductivity is read out^° from Eqs. fHT]) and P2l) : 

/2(ik df 
{2-nf g^ ^[-^(^^)]^M(k)^^(k), ^i,v = x,y,z (43) 

Using Eq. (l2^ and the definition of velocity operator = ^gf-, Eq. (l40|) is changed into 



<'M = S E "w..„.EE(i-«-)(^^4x^+ ' ^ ' 



Im{[ / rf^xxXX/Ji^ / d^^ixlv^Xd-XdV^^xt)]} (44) 
Jn ^ Jn 

If we interpret as the energy-dependent relaxation time t{Ei^) caused by inelastic 

scattering of phonons (in a given MD step, scattering is caused by deviation from crystal; 
to refiect various vibration states and electron-phonon scattering, averaging over many MD 
steps is necessary^"^) , n/^(l — nd){ ^^ -E^+huj + Ei -Ei-h^ ) Ie' ^^"^ states as integral 
over Brillouin zone: ^ T^hi^-.i^^ ^hh-iN, J2a=i ^ / (only in crystal, one can 

use k-points in the reciprocal space to characterize states), where ni^ is removed from W. 
At a; = 0, as expected, Eq. (j44|) is reduced to the semi-classical result (j^3l) . 
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The sinut terms terms in Eq. (l39|l lag 90" behind the phase of external field. The contri- 
bution to the imaginary part of conductivity is 



2fc ^ 

(2). s e h 



1 . % ■[ 2 



MiXdlMxiJ [ d'xix,^ - Xl^)] (45) 

Jn ox a ox a 

it is interesting to notice that o"^^^(0) = 0. 

To first order in the field, the second term in Eq.([B]) is 

- -A(r, i)^{Wi...w^}(r, t) = ^^{w,...w^}(r) sina;t (46) 
From Eq.Q, to 1st order of field, the 3rd term in Eq.Q is 

J]^sin.;t5(r-W„) (47) 



a 



Using Eq.([H]), the imaginary part of conductivity of the electrons + nuclei is 

„2„e n'^n^ 

.i.M-SM.^.E^ (48) 



Mp, Qp and rip are mass, effective charge and the number density of the p*^ species of nuclei. 
The last two terms are contributions from free charges^. 

We cannot use Eq. (|371) in two situations: (1) degenerate states with low frequency field, 
= and — 0; (2) external field and two groups of levels in resonance: Udc — c<j = or 
Udc + uj = 0. In these situations, Eq. fl36|) leads to aa(t) ~ t. 

2. Interaction of degenerate states with a very low frequency external field 

For a group of M degenerate states {xd,y, cr = 1, 2, ■ ■ ■ , M), the mutual coupling is much 
stronger than the coupling between one member and the states with different energy. The 
general evolution equation 

^n^e-^^*/^ = ^afc\/,,e-^^'=*/^ V,, = [ drx*HjUt)Xk (49) 

7., 



is simplified to 

_daj{t) 



ih 

dt 

k 



^ak{t)Gjk COS ut, Gjk= drXj{-eEo-r)xk j,k = di,d2,- ■ ■ ,dM (50) 

7^ ^ 
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Introduce new variable s = sin tut, Eq. (l50|) becomes 

thw^^ = J2^kis)G,,, (51) 

k 

Notice Gjk does not depend on time, taking Fourier transform 

aj{s) = j dpapje^'P', j = di, c/a, ■ ■ ■ , c?a/ (52) 
in both sides of Eg. (1511) leads to 

apj{Gjk - h^jp5jk) = (53) 

k 

Because (Gjk) is Hermitian, its eigenvalues pi,P2, ■ ■ ■ ,Pm are real and the matrix (a^^^^) is 
unitary. The M special solutions of Eq. fISUl) are 



M 



40)(t) = e-^'^^^-'^*^a,^,„X.., /x = l,2,---,M (54) 

a=l 

The matrix elements of perturbation (— eEo-r) relative to the new zero order wave functions 
are diagonal: 

j drxJ°)*(-eEo ■ r)x?? = Gp^^J,,, /x, z/ = 1, 2, ■ ■ ■ , M (55) 
Because matrix (op^^^) is unitary, one has 

M M 

EW? WVx?:)*(t) - x'iT{t)VxZ\t)] - Vxl - xlVXdJ = (56) 

(1=1 a=l 

From Eq. (l3T]) . the zeroth order correction to degenerate states does not contribute to con- 
ductivity. 

Unlike the KGF, where dc conductivity is obtained either by extrapolating from optical 
conductivity^ or by writing a separate code for zero frequency^i^'iS, the present ac expression 
includes dc expression in the obvious way. One may notice when u = 0, = and j^2 — 
(cf. Appendix A). Excepting the factor (1 — n) which would not appear from single-electron 
stationary perturbation theory, when u; — > 0, is reduced into the second sum of Eq. fl30l) . 
(jc2 + Jo) is reduced into the third sum of Eq. (l30|) (cf. Appendix A). 

In Appendix A, we show that the contributions to current density or conductivity from 
the degenerate states is finite. In the single-electron states, the number of degenerate states 
in each degenerate manifold is much smaller than the total number of non-degenerate states. 
The conductivity from the coupling among degenerate states can therefore be neglected. One 
can use Eqs.( l40|) and (H5|) . in which only counts the coupling among non-degenerate states 
and the coupling between degenerate states and non-degenerate states. 
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3. Resonance with external field 



Suppose in {xh iXhi' ' ' i XIn^ ) there is a M— fold degenerate manifold 
(Xmi,Xm2, ■ ■ ■ ,XmM) ^nd a M'-fold degenerate manifold (x^, Xna, ■ ■ ■ , Xn^v^,) which 
are nearly in resonance with a finite frequency u external field: e'^ — En^ = h{uj + e), 
e<<co'ore = 0. The coupling with other non-resonant states can be neglected. If we only 
consider the interaction with the smallest oscillating frequency, 

the general evolution equation 
da 

^^-^ = 5^V;^,n,(t)a„„ i = l,2,---,M; fc = l,2,---,M' (58) 

k=l 

is simplified to 

'f^^ = J2^rn,n,e'''a^, (59) 

k=l 

Similarly if one only takes the interaction with smallest oscillating frequency 

V (f) ^ F* p«*('^"fe™j+'^) — p* p-it^ (ar)) 

the general evolution equation 

da 

^^-^ = i^)"'rn, (61) 

is simplified to 

, M 

^^-^ =Y.^kn.--'''-r., (62) 

Introduce new functions hn^,'- an^ = feufcC"**" {k = 1,2, ■ ■ ■ , M'), Eq. (l59|) becomes 

da ^' 

^^^f = Y F^,n,bn, (63) 

k=l 

Eq. (!62!) becomes 

M 

ihbn, = -hebn, + F^^n^am^ (64) 



Eqs.(ES]) and can be rearranged into 



ih^V = RV (65) 
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where 

R=( ^ I (66) 

is a (M + M') X (M + M') matrix, is zero matrix, I is the unit matrix, elements of matrix 
B are given by 

Bjk = Fm,n„ J = 1,2,---,M; A;= 1,2,--- ,M' (67) 
is a (M + M')— column vector, its transpose is 

We are looking for special solutions of Eq. (jHSD in the form: 

<,W=<e^*"% J = 1,2,---,M; 6^^(t) = 6f^e^*"^ A: = 1, 2, ■ ■ ■ , M' (69) 
The column vector Vq with 

(^^,^.anspo.e ^ ^^^0^^ ^^0^^ . . . ^ ; < , 6f„ ■ ■ ■ , 6^^, ) (70) 

is the eigenvector of R belonging to eigenvalue hag. Then (M + M') special solutions of the 
time-dependent single-electron Schro dinger equation are 

M M' 

xT\t) = J2 <e^°^*Xn.,e-^^™*/'^ + ^ <e^("'-^)*X.,e-^^"*/^ g = 1, 2, ■ ■ ■ , M + M' (71) 
j=i k=i 

and the general solution can be obtained from linear combinations. Because R 
is Hermitian, its eigenvalues {hag, q = 1, 2, ■ ■ ■ , M + M') are real, matrix C = 
(«mV ■ ■ ■ . «mM5 ^ni' ^n2' ' ' ' ' ) Unitary (g is index of row). 
If we use Eq.f l7ip and notice that C is unitary, one has 

M+M' M M' 

9=1 j=l k=l 

(72) 



From Eq. lUSl) . the contribution to current from two groups of states in resonance with an 
external field is zero. The artificial poles in the case of resonance in Eq.( l37I) are removed. 
For a mechanical oscillator, if we input energy in a resonant way and do not take out 
energy, the amplitude of the oscillator will increase indefinitely. The situation for two 
groups of resonant levels is different; the system absorbs the external field while stimulated 
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emission also occurs. The material and field are in absorption-emission equilibrium, so that 
no singularity occurs. The two 5 functions (they originate from first order correction of 
wave function) in Greenwood formula come from the long time limit, and are not caused by 
resonance. In Appendix B, we show that the contribution from resonant states is finite. 



C. Comparison with Greenwood formula 



Both the present work and the Greenwood derivation require the use of perturbation 
theory: 

ad{t) = -^Fdc [ ' dt'[e'*'("''^-"^ + e'*'(-''^+-)] « 1 (73) 
Jti 

This means that the interaction time r = (^2 ^ ^i) cannot be too long: 



r = (t2-ti)«| (74) 

The Greenwood derivation also requires that the transition probability per unit time be 
defined: 

lim — — — = n6( ) and lim — , , , = nSi ) (75) 

That is, the interaction time r should be long enough 

2 2 
T » and r >> (76) 

to allow the two limits in Eq. flTSl) to be taken. The law of conservation energy (of field + 
matter) can be verified by means of two measurements only to an accuracy of the order of 
h/At, where At is the time interval between the measurements^, i.e. the interaction time 
r between field and matter. Since in present work we do not need probability per unit time 
(i.e. long time limit), the energy conservation delta function will not appear. 

For a large system with continuous energy spectrum, Eq.( |76l) contradicts Eq. (!74l) for 
close levels when u 0. Since a transition with small Udc makes a large contribution to 
conductivity, the dc conductivity obtained from the KGF is problematic. The derivation in 
this work does not need transition probability, therefore does not need condition (!76l) . and 
is self-consistent. Numerically Eq.(l40l) and Eq. (l48l) get rid of the delta function, and do 
not require an artificial broadening. 
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III. ROLE OF MANY-ELECTRON STATISTICS 



In this section, method (2), many-body perturbation theory, is used to compute the 
conductivity. We take a static field as example and apply to concrete examples. 

A. Intrinsic semiconductor 

1. K- electron excited state 

Label the states in valence band from low energy to high energy as v^^, ■ ■ ■ ,V2,Vi, the 
states in the conduction band from low energy to high energy as Ci,C2, - ■ ■ ,cjv^- For an 
intrinsic semiconductor at T —OK, the valence band is full, and the number of states is the 
number of electrons. The system is in its ground state 

At T > 0, various excited states appear. First consider one electron pumped into the 
conduction band from the valence band. A one-electron excited state ^vjCk is constructed 
from Aq by replacing Vj with c^. There are manners in choosing Vj. There are A^e manners 
in choosing Cfe. The single-electron energy spectrum is very dense, there are many combi- 
nations of the choices of initial valence state and the final conduction state, the degeneracy 
among 1-electron excited states is high. 

Since any observable, including current, are bilinear forms of a A^g-electron wave function, 
the order of states in the Slater determinant A^^c^ does not matter provided that we 
maintain the same order in A*^^^. The probability of state ^.^^ck relative to the ground state 
Ao is 

f/.^e, = exp{-[(E„, - E,^) + E, + - E,^)}/{kBT)} (78) 

where Eg is the band gap. 

If two electrons are pumped from valence band to conduction band, a 2-electron excited 
is obtained from Aq by replacing {vj^,Vj^) with (cp^jCp^). The relative 
probability of state ^vj^vj^cp-^cp^ to ground state Aq is 

Uv,,v,,cp,cp, ^ eM-KEvi - E,.^) + {E,, -E,.^) + 2Eg + {E,^^-E,J + {E,^^-E,J]/{kBT)] 

(79) 
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In Ao if K electrons are excited to the conduction band from valence band: state 
Vj^ is replaced by state Cp^, state Vj^ is replaced by Cpj, ••• , state Vj^^ is replaced by 
state Cp^, a K-electron excited state ^Vj^vj^-vj^^^cp-^cp^-cpj^ is obtained. The probability of 
V -v c c -c relative to An is 

K K 

Uv,,v,,-v,^cp,cp,-cp, = exp{-[^(E,, - E,J + KE, + Y.^E,.. ' ^cJ]/(^sT)} (80) 

a=l a=l 

The absolute probability W^-c^ of 1-electron excited state I^vjCk is 

jk ^3\^32"''"iK^'^Pi-'^P2"''^PK 

It is easy to verify for low T that exp(^^^^^) << 1 and exp( '^°~^^ ) >> 1 (they are satisfied 
even in several thousand K), so that one has 

W,^,, = [1 - f{E,^)]f{E,,) (82) 

where 

are Fermi distribution functions of valence states and conduction states. 

The absolute probability ^vj-^vj^cp^cp^ of 2-electron excited state ^vj^vj^cp^cp^ can be ob- 
tained similarly 

^vi^vj^cp^cp^ = Uv^,vj^cp,cpjZ = [1 - f{E^^^)][l - f{E^^^)]f{E^^)f{E^^) (84) 
The absolute probability ^vj^vj^-vj^^cp^cp^-cpj^ of K-electron excited state 

^Vj^vj^-Vjj^;cp^cp2-Cpj^ is 

K 

^Vj^Vj^-Vjj^Cp^Cp^-Cp^ — Uvj-^Vj^-Vj^Cp^Cp^-Cp^l Z = — f{Evj^)]f{Ecp^) (85) 

a=l 

2. Zero dc conductivity at T=OK 
Because the interaction with a static field 

Hint ^~Y.eE-Vm (86) 



m=l 
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is a single-particle operator (separable for coordinate of each particle), the ground state only 
couples with 1-electron excited states 

(A„^,,J - J] eE ■ r^lAo) = / (iri4(ri)(-eE ■ ri)t;,(ri) (87) 

m=l 

The change in ground state Aq by external field only includes 1-electron excited states 
^,(,) _ ^ (A,,,J - E^^L, eE ■ r^lAo) ^ _^ {c,{l)\e^ ■ v,\v,{l)) 

(88) 

If we take Eq.([6]) and effect the multiple integral, the current density is 



(89) 

^ (cfc(l)|eE ■ ri|fj-(l)) f J t t \V7 *t \ *t \V7 I \\^ 

+ 5 (gc. - + (g.. - + I, *("=('^)^-"^ - (^)^-''(=>» 

Using Eg. (1221) . one can read off conductivity. Because the external field is much weaker than 
the atomic field, the numerator is much smaller than the energy gap Eg (this will become 
more obvious in next section), the change Aq^'' in wave function Aq can be neglected, and the 
dc conductivity is negligible at T=0 in an intrinsic semiconductor. The coupling between 0- 
electron to 1-electron excited states can be viewed as a cross band transition, its probability 
is not exactly zero, but is extremely small. One may neglect the existence of conduction 
band: electron cannot be accelerated when valence band is full. To accelerate an electron 
in ground state, one has to go from valence band to conduction band. The probability is 
negligible for an external field which is much weaker than atomic field. 



3. Conduction from one- electron excited states 

Because Hint is single-particle operator, a 1-electron excited state could couple with 
ground state, 1-electron excited states and 2-electron excited states. The energy difference 

between a 1-electron excited state and a 2-electron excited state is at least energy gap Eg. 
The contribution to current density from this coupling is small. So does the coupling between 
a 1-electron excited state and Aq. 
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Since Hint is a single-particle operator, there are only two types of coupling between two 
different 1-electron excited states: A^^.cfc ^ ^v^cy or A^.c^. ^ A^^/Ck- The energy difference 
between such pairs of 1-electron excited states can be small if states c^/ and Ck {vji and Vj) 
are properly chosen. Their contribution will be much larger than the coupling between a 
K-electron excited state and a {K ± l)-electron excited state. 

The change in A^.c^. caused by a static field is 

a'(i) _ (cfc/(l)|eE-ri|cfc(l)) ^ (ty (l)|eE ■ ri|i;j(l)) 
^^^-^ - ^ , - E^^ + ^ E^,-E^. ^^^"^^ 

By appealing to Eqs. f[T^ . fl22]) and ([23]), the expression of conductivity in momentum 
representation is 



jk k 
d 



doc r\ f 



[Ck'[r)— Cfc^rj 



The accelerated hole in the valence band and the accelerated electron in conduction band 
contribute most to the conduction, the coupling between i^-electron and {K ± l)-electron 
excited states contributes much less. We have proven that degenerate states act like non- 
degenerate states, cf. Eq. (l30!) . All the denominators in Eg. (191]) are not zero. Except 
delta functions, Eq. flOT]) . the contribution from 1-electron excited states, corresponds to the 
ordinary Greenwood formula. 

In the standard procedure applying KGF^ 

dEa{E)[-M, a{E) = ^ J] | (nb,.|m) - E)6{E„, - E) (92) 

one broadens delta function by a Gaussian 

,(S„_B)«!iPhW.Z^£)!M (93) 

Av 277 



Numerically, this procedure is equivalent to replace whole series about {E^., — E^.)'"^ in 
Eq. (19T]) with several large terms, each of order of A"^. There are two relevant energy scales 
in the problem: fc^T, and characteristic energy level splittings near Ep. The choices of 
A is thus somewhat arbitrary. On one hand A should be order of ksT to reflect thermal 
environment. However /c^T is a too small choice of A for room temperature but may be 
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too large for a high temperature. On the other hand, A should be order of or smaller than 
the eigenvalue splittings near Ep. This choice depends on the size of a structural model and 
also depends on how many k-points one wishes to use. Thus KGF depends on a fortunate 
choice of A, or requires some other extrapolation scheme to u; = 0. Eg. (191 1) or Eq. fHU]) does 
not suffer from this problem. 



4- Conduction from 2-electron excited states 

Although a 2-electron excited state may couple with a 1-electron excited state or a 3- 
electron excited state, the energy differences are at least energy gap Eg. Later, we only 
consider the matrix elements between two 2-electron excited states. The 1*^* order correction 



to -/Vti ■ 7t - n IS 

^31 ^32 PI P2 



^ (cp/(l)|eE-ri|cp,(l)) ^ (cp^(l)|eE ■ ri|cp,(l)) 



where we keep the single-electron wave functions in each Slater determinant in a fixed order. 
Substitute Eq. (!M|) into Eq.([6]) and effect the multiple integral, the macroscopic current 
density is then: 

ii J2(>ii) pi,P2{>pi) 



^ (i;,-'(l)|eE- ri|t;,-,(l))* f 



{v^,^{l)\eE-v,\v^,{l)y 
^ (t;,-'(l)|eE-ri|f,-Jl))* f 



Ev,, - 




(^i^(l)leE- 


ri|^i2(l))* 


^^.^ " 


- E^, . 


(c,,(l)|eE- 


ri|cp,(l))* 


^S' - 

p\ 


- Ec^^ 


(cp^(l)leE- 


ri|cp2(l))* 


^v, - 


-^c^2 


(t;,.(l)|eE. 


■ril^ji(l)) 



^ (c„'(l)|eE-ri|cp2(l))* /■ 
^ (t;j/(l)|eE-ri|t;,,(l)) /" 
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Pi(7^Pl) 



(t;,.(l)|eE-ri|t;,,(l)) 



P2(^P2) 



^-.^ - 






(cp;(l)|eE- 


ri|cpi 


(1)) 








(V,(l)|eE- 




(1)) 









dr(c,,(r)VrC;^(r)-c;^(r)VrCp.(r)) 

dr(c,, (r)VrC* (r) - c* (r)VrC,.,(r))} (95) 



Using Eq. (l22|) . one can pick off the conductivity from coupling between 2-electron excited 
states. 



B. Metals 

In a metal, the conduction band is partly filled. Relative to the Fermi surface, holes and 
electrons are in the same conduction band. The energy difference between hole and electron 
always can be taken as small. Beside the non-existent energy gap, the difference between a 
metal and a semiconductor is that the number of carriers ~ ^"-^e in the former is greatly 
larger than that in the later. It is easy to check 



/ 



dridr2dr2 



l[{2) l[{3) 
I'.il) l'^{2) /^(3) 

^3(1) ^3(2) m 



(ri + r2 + rz)- 



h{l) h{2) h{3) 

kil) km his) 

hil) h{2) kiS) 



(96) 



there are 9 terms, each term has 2 sub-terms (they form a determinant). By induction 
method, one finds 



m=l 



j,k=l 



where Dj^ 



is a {Nf. — 1) X [Nf. — 1) determinant, each element of D is a Kronecker 
delta symbol. The row indices are (Z1/2 ■ ■ - ^n^) in which Ik is removed, the column indices 
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are {l[l2 ■ ■ - I'nJ iii which /'• is removed. For example 



D 



hi' 



■V I' ■■■/' 







^lA ■ 


'I'jVe 






^hi'^ ■ 






^lA 


^lA 






^InA2 


^InAa ■ 


■ ^InAn 



(98) 



The first order change in the iVe— electron wave function is 



A, , , 

Next, substitute Eq. 
tivity: 



jk 



A/' /' .../' Dnh 

'l'2 'iVe J 



(99) 



El, - El, 

into the expression of current density, one finds the dc conduc 



'l'2 'JV, 



a 



rriQ 



\j+k 



hh-'-lNe 



p,q=l 

(/'j^i2'"'jV„''l'2---'iVe) 



(''l'2'"'jVe''l'2---'iVe) 
pi} 



jk 

dx,, 



where Dpq ^ ^ ^ is a (A^^e — 1) x [N^. — 1) determinant, each element of which is a 
Kronecker delta symbol. The row indices are (/'^/2 ' ' ' ^jVe) which is removed. The column 
indices are (Z1Z2 ■ ■ - ^n^ in which Ip is removed. 

1 



W^z.Z2- 



Ei.-Efs 

l + exp(^^. 



;ioi) 



is the appearing probability of Ai'e— electron state A^^^j...^^^ 
distribution, it is easy to check 
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Using the definition of Fermi 
E-Ef (102) 



g/3e _^ I 

introducing an electron in a state above Fermi surface {E > Ep) is equivalent to introduce 
a hole below Fermi surface {E < Ep). It is clear from Eq. OlOll) . the states around Fermi 
surface contribute most to conductivity, as expected. 

In a semiconductor, due to the energy gap Eg, the appearing probability of a K-electron 
excited state includes a factor ^-^Eg/ksT^ calculate the conductivity, it is enough to 
restrict attention to the excited states with few electrons. In a metal, the conduction band 
is half-filled, and there exist many low-energy excited states. One must count all electrons 
although only a shell close to Fermi surface makes an important contribution. 
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C. Homogeneous doped semiconductor 

In a weakly n-type doped semiconductor, there are substitutional atoms or interstitial 
atoms. The energy levels of the former lie just below the bottom of conduction band; the 
energy levels of the later lie above the Fermi level /Xj of the intrinsic matrix. There are three 
contributions to conductivity: (1) electrons from substitutional donors, given by Eq. fllOOl) . 
EiT' is the chemical potential of doped material in Eq. fllOip : (2) carriers from interstitial 
atoms; (3) carriers from intrinsic matrix. Parts (2) and (3) can be calculated by Eq. (!9T!) . 
Similar consideration is applicable to weakly p-type doped semiconductors. 

IV. SUMMARY 

We discussed some foundational issues with respect to computing the conductivity, and 
improved the Kubo-Greenwood formula by computing dc and ac conductivity from current 
density, in which only the amplitude of probability, and not the transition probability itself 
is used. Eqs. (1401) and (jUj) are key new contributions of this paper. In this method, the 
expression of dc conductivity is extracted from the ac conductivity in a direct way. We found 
that (1) the contribution from the states which are near or in resonance with finite frequency 
external field is finite; (2) the contribution from degenerate states in low frequency or zero 
frequency external field is finite; (3) the energy conserving 5-function does not appear in the 
improved expression, thus one can avoid artificial numerical broadening. In the formulation 
of many-body perturbation theory, i.e. "method (2)" for calculating the current density, the 
many-electron statistics is displayed explicitly. One example is that the dc conductivity of an 
intrinsic semiconductor at T=OK is zero. For an intrinsic semiconductor, Kubo-Greenwood 
formula is the contribution from 1-electron excited states. 

V. APPENDICES 

A. Conductivity from degenerate states in a low frequency external field 

We first compute the first order correction to Xp'^J{t): 
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Now the zeroth order wave functions are 

... y, ... ■ y'(0) yXO) . . . y'iO) 

At an initial moment, one electron is in state Xp'^J: ap^(— oo) = 1 and other coefficients 
are zero. If the interaction time with field is not too long, ap^{t) is dominant. For a non- 
degenerate state Xjj ^^Pf!ji^) determined by 



daj{t) ^ (J'-Jt I f,-i'-Jt\n (-i'\J'^^r,„t 

dt ~ 2n '"^^ + ^ "^^ 



-GjpAe''^' + e-'"*)ap rt)e'"^''M* (104) 



The solution which satisfies initial condition aj(— oo) = is 

= + -], 5^0+ (105) 

For a member of the degenerate states Xp^,, (/w' 7^ A*); '^p^iif) satisfies 

.^^^e-'^/^ = E a.(t)G,^,,e-'^^^/^ ^" (106) 

k 

and initial condition a^^, (— oo) = 0. Index k in RHS of Eq. fll06p runs over non-degenerate 
states only. Because Eq. fl55|) . no coupling among {xpm''} in Eg. (11061) . cip}p^i{t) is given by 



a. 



4;^2 Z-j m' ^^^^^ + ^ ^^^^^^ + 2u; + 2(5 t^^p,, + ujp^,p^ + id 



+ ^ ^ + ^ ^} (107) 

cjfcp^ - uj ujp^,p^ + lb ujkp^ - uj ujp^,p^ -2uj + id 

When — > 0, all the denominators of Eqs. (11051) and (11071) are non-zero. For degenerate 
states in a zero frequency external field, the artificial singularity of perturbation formula 
(137|) is removed. 

Combining Eqs. (ll03p . (I105p . (11070 and fl^ . the macroscopic current density can be 
written as 

f (r',t) = jLn(r',t) + j:(r',t) + j:(r',t) + U^\t) + U^'^t) + jg(r') (108) 

where j^onl^'? ^^^^ contribution from non-degenerate states, and is obtained by replacing 
Yl!l=i by (sum over only non-degenerate states) in Eq.( l39il . 

M 

j:(r',t)=cos^t^-^ W^M.-^..EE(1-^^) (109) 
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is the component with cos ujt factor. 

M 

j:(r',t) = sina;t-^ J] W^W.-/.. E " ^^■) (^O) 



r' 



1 



^] Re G,,^ / ds{x,V.4T - 4r V. 



is the component with sin cut factor, 'fc fs come from couphng non-degenerate states 
with the degenerate manifold. 

M 



j:,(r',t) = cos2o;t--^ E ^^^.-^e E " V) E (^^l) 



is the component with cos 2ujt factor. 

M 

■ \^ w, , , 



rAr\t) = sm2ujt-^ E ^M.-UeEE(i-v)E (112) 



11 11 



ReGp^,.G,p^ I ds{4y.4^>-4^>v.4^,) 



^kp^ - ^ ^p^>p^ - cjfcp,, + uj ujp^,p^ + 2uj 
is the component with sin2co't factor. 

M 1111 

jo(r') = T^ E w^M....eEE(i-ME(^T — — + ^ —) 

(113) 

Im Gp^,,G,p^ I dsi4:^*v.xTJ - 4')^^4T)\ 

is the component without time variation factor. Using Eq. fl22|) . one can easily identify the 
conductivity. 

In Eqs. fll09p - flll3p . the summation over j or k is not restricted to (/1/2 ■ ■ ■ ^atJ: it extends 
to all single-electron states. js2 come from indirect coupling among {xpm'*' ~ 

1, 2, ■ ■ ■ , M} through non-degenerate states. Interaction with external field appear twice 
in Eg. (11071) . three new time factors cos2c(jt, sin2c(jt and 1, which are different from the 
original time factors e"**"^ and e**"^ of the external field, come from the 2^^ order harmonic 
generations. 
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B. Conductivity from resonant states 



To compute the contribution of two groups of resonant states to the current, we need the 
I''* order wave function 

xT\t) = 5^af )e-^^*/\. + $^<i^)e-^^"*/\„, + J2<^"^''''^'Xm, (114) 

s k j 

where s indicates the states which do not belong to the upper and lower degenerate 
groups, j scans over the upper group {Xmi,Xm2^ ' ' ' jXtom)? ^ scans over the lower group 
iXni, Xn-ii ■ ■ ■ , Xnj,.j,)- Supposc initially the system is in the g*^ mode of the resonance states, 
the first order evolution equation is then: 

j 

+ Et^-^^^*^""^"^^ + ^4.e^*^-'^^+-)]a(°) (115) 

k 

The solution of Eq. ( 1115^ is 



a '—^ Usrn: — UJ -\- ttq — e ^srik + + CHg — 6 

For a member of the upper group, the first order probability amplitude is determined by: 

oo^ ii(2c.+.) (0) I V e"**"^ + F* e^*^la(°) fll?) 

+ Et^-.^^'*^"'"'''""^ + i";n^,e'*^"'"^■^''"Vi'^ J = 1, 2, ■ ■ ■ , M 
Using Eq.f lll6p . the solution of eq.f lll7p is 

fe=i 



qO 



^f^2l^ ™.-tZ^^ _ ^ + _ 2c.) ^ "^^-^ (c. + + «,)«,^^'"/ 

„i{aq—u))t A{aq+u))t 

+ S^\F - + F* - 16'?° I 
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1 ..-^ P«*"9 pi{aq+2uj)t 

S J J J 

For a member of the lower group, the first order probabihty amphtude is determined by: 

^«a^ - V F p-it(2a;+e) (0) , F p p-»t<^ , P^^'^lflW 

+ Xlt^"*^^'*^''"''"^^ + ^.nfee'*^''"''^'^^]«i^^ = 1, 2, ■ ■ ■ , M' (119) 

s 

Using Eq.dlin]), the solution of Eq. (|TT9l) is 

i=i ^ 

(120) 

+_Vf iVfF - +F* - la'^° 

^K^l-^ "'^"^^^ ^™^(^^^^, _^ + a,)(a,-3^-e)^ + 

^ pit{aq—2ui—e) pit(aq—t) 

^ {uJsny - UJ + aq - (:){aq - 2UJ - e) ^ {Usny + O; + - c) (ftg - e)^^"'='^ 

1 pit{aq—u:—t) it{aq+u)—t) 

iX^fF - h F* - 

^ s j [UJsm, - + <yq>{<yq - UJ - e) ' (c^^m, + + ("g + - c) ' 
^ gii(aq-e) pit{aq+2ui—e) 



+ - e) (ftg - e) {usny + uj + ag - e){aq + 2uj - e) 



Substituting Eqs.([Il6l), ffTTSD and ffM into Eqs. ffTTil) . by means of Eq.(l3l]), one can 
find current density. The full formula is too long to write out, we only write down the 
contribution from coupling between non- resonant states and resonant states: in Eq. flll4p 
only keep the first term. The current density with time factor cos cut is 

jc(r,t) = ^cosa;t /" ds ^hh-iN^Y.^^-''^^ (^21) 



Im{ J]a^%^°;(x™, Vx: - xyXm^m + F 



1 

IS 



The current density with time factor sin tut is 
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jsir,t) = -^smujt f ds J2 (122) 
Ini[J]a^O;a^° (x.Vx™^. + X™, Vxs)[F.m,, ^— F*^ ^ ] 

^ ' ' ' ' ' ^sra,, - ^ ^ OLq ^ UJsm,, + + 

The current density with time factor cos 2ujt is 



"'"''tUsn. + + a„ - e 



1 



The current density with time factor sin 2ujt is 



js2(r,t) = -^sin2a;t /" rfs W^w^-^^, $^(1 " ^s) (124) 

imK%f;(x„., vx: + x:vxrn,)F„,, 

+<C(x:.vx. + xsVx:jF„ 



1 

' + + - e 

1 



t^smj - W + 



The current density without time variation factor is 



1 



ReK° 6t(x™, vx: + x:vx^,)f;„, — ^ 

. — i 1 



+<C(xn.vx:+x:vxnjF, 



The contribution to current from two groups of resonant states is finite. Using Eq. (l22|) . one 
can again read off conductivity. 
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